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An annihilation signal of dark matter is searched for from the central region of the Milky Way. 

Data acquired in dedicated ON/OFF observations of the Galactic center region with Fl.E.S.S. are 
analyzed for this purpose. No significant signal is found in a total of ~ 9 h of ON/OFF observations. 

Upper limits on the velocity averaged cross section, (av), for the annihilation of dark matter particles 
with masses in the range of ~ 300 GeV to ~ 10 TeV are derived. In contrast to previous constraints 
derived from observations of the Galactic center region, the constraints that are derived here apply 
also under the assumption of a central core of constant dark matter density around the center of the 
Galaxy. Values of (av) that are larger than 3 • 1CV 24 cm 3 /s are excluded for dark matter particles 
with masses between ~ 1 and ~ 4 TeV at 95% CL if the radius of the central dark matter density 
core does not exceed 500 pc. This is the strongest constraint that is derived on (av) for annihilating 
TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution 
in the search region. 


INTRODUCTION 


The formation of the large scale structure of the uni¬ 
verse as well as the dynamics of galaxy clusters and in¬ 


dividual galaxies strongly suggest the presence of dark 
matter on the respective length scale [l|. Many exten¬ 
sions of the standard model of particle physics predict 
a stable particle without electromagnetic coupling whose 
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presence can account for the missing mass that is appar¬ 
ent in astrophysical environments [l|. The annihilation 
of dark matter particles is expected to produce photons 
with energies up to the mass of the dark matter particles 
0. The detection of 7 -rays from a given direction can 
thus indirectly probe the presence of dark matter parti¬ 
cles along the corresponding line of sight. 

The central region of the Milky Way is of particular in¬ 
terest for indirect searches for annihilating dark matter 
because the squared dark matter density integrated over 
the line of sight towards the target region (i.e. the as¬ 
trophysical or J-factor) is expected to be large 0 . The 
J-factor for observations of the Galactic center region de¬ 
pends strongly on the dark matter density distribution 
within the Milky Way. Simulations of the dynamics of 
the dark matter content of galaxies predict to universal 
dark matter density distributions. Towards the center 
of the galaxies, the influence of baryons on the distribu¬ 
tion of dark matter is not yet resolved. The formation 
of pronounced density cusps towards the center of galax¬ 
ies (0, jl5j]) and, more recently, the prediction that 
the dark matter density in the central few hundred pc is 
almost constant @, [)|, [fij] have been discussed. The lat¬ 
ter prediction of an almost constant dark matter density 
in the central region of the Galaxy is considered in this 
work. 

The current strongest constraints on the velocity aver¬ 
aged cross section for the self-annihilation of dark mat¬ 
ter particles with masses in the range of ~ 400 GeV to 
~ 10 TeV result from a search for an extended emission 
of 7 -rays in the central region of the Milky Way with 
H.E.S.S. [7j. However, the constraints apply only if the 
dark matter density distribution in the central ~ 500 pc 
of the Milky Way is cusped (see dotted lines in Fig. [TJ. 
An alternative search for the annihilation of dark matter 
particles is presented in this paper and strong constraints 
on (av) are derived without the assumption of a dark 
matter density profile that is cusped in the central 500 
pc of the Milky Way. 

THE HIGH ENERGY STEREOSCOPIC SYSTEM 



FIG. 1. Dark matter density as a function of the distance 
to the Galactic center. The parametrizations for the cusped 
Einasto and NFW profiles are taken from iijj]. The cored 
Einasto and NFW density profiles follow the respective cusped 
profiles at distances to the Galactic center that are larger than 
the core radius of 500 pc. 

constructed. The construction of the background region 
must be performed such that the ratio of the acceptance 
for background events in the signal and the background 
region is known from instrumental characteristics. The 
definition of the background region enables a compari¬ 
son between the number of events that are detected in 
the signal region and the number of background events 
that are expected in the signal region. The acceptance 
for background events of H.E.S.S. is in general strongly 
influenced by atmospheric conditions, the pointing zenith 
angle and the night sky background in the observed field 
of view. See [9| for more information on the H.E.S.S. 
experiment. 


ON/OFF OBSERVATIONS OF THE GALACTIC 
CENTER REGION WITH H.E.S.S. 


The High Energy Stereoscopic System (H.E.S.S.) is 
an array of imaging atmospheric Cherenkov telescopes 
(IACT) in the Namibian Khomas Highland. IACTs 
detect the Cherenkov light emitted by electromagnetic 
showers that are induced when primary 7 -rays interact 
with air nuclei in the Earth’s atmosphere. Charged cos¬ 
mic rays also induce showers in the Earth’s atmosphere 
and constitute background for the IACT detection of 
7 -rays. Cosmic ray background events that cannot be 
suppressed during the analysis of H.E.S.S. data (see i) 
are typically treated with a background subtraction tech¬ 
nique 13]. The background subtraction relies on the def¬ 


inition of a signal region for which a background region is 


The ON/OFF observation mode (see also 0, 0) 
refers in this study to a special observation strategy 
where a background (OFF1) region, the signal (ON) re¬ 
gion and another background (OFF 2 ) region are observed 
consecutively for 33 min each. Figure [2] shows the ob¬ 
served regions in galactic coordinates. The signal region 
has a radius of 2° and centers at l = 1°, b = —0.7° 
in galactic coordinates or a = 267.7°, S = —28.4° in 
equatorial coordinates (J2000). The centers of the two 
background regions have a symmetric offset of ±35 min 
in right ascension to the signal region center. The two 
minute difference between the right ascension offset be¬ 
tween the signal and background regions and the obser- 
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FIG. 2. The signal region close to the Galactic center and 
the two background regions with a symmetric right ascension 
offset of ±35 min to the signal region. A right ascension - 
declination coordinate grid is overlaid. The exclusion of the 
Galactic plane (|6| < 0.3°) in the signal region and the regions 
with a ±35 min offset in right ascension to |6| < 0.3° in the 
background regions are visible. 


FIG. 3. Line of sight integral over the squared dark mat¬ 
ter density as a function of the angular distance between the 
direction of the line of sight and the direction of the Galac¬ 
tic center. A 500 pc core radius is assumed for the cored 
dark matter density profiles. Overlaid is the range of angu¬ 
lar distances to the Galactic center covered by the signal and 
background regions of the OFF1/ON/OFF2 datasets. 


vation length allows for a transition time between the 
observations. The ON/OFF observation pattern allows 
the equalization of the azimuth and zenith angles that 
are covered by array pointings in each of the obser¬ 
vations. Differences in the acceptance for background 
events which result from differences in the zenith angle 
array pointing range can thus be neglected. The time dif¬ 
ference of 35 min between the observations is a compro¬ 
mise between the demand for small atmospheric changes 
(i.e. small time differences) and a large offset in right 
ascension (i.e. large time differences). Two background 
regions are observed, to better control residual imbal¬ 
ances in the acceptance for background events between 
the observations. 

Figure [3] shows the J-factor for a given line of sight as 
a function of the angular distance, 9 , between the direc¬ 
tions of the line of sight and the Galactic center. The 
J-factor is proportional to the expected number of dark 
matter annihilation events in the respective direction. 
The 6 angle ranges that are covered by the signal and 
background regions in the OFF1/ON/OFF2 observations 
are indicated in Fig. [3] It is concluded from this figure 
that the expected number of dark matter annihilation 
events is larger in the signal than in the background re¬ 
gions when the radius of the core of constant dark mat¬ 
ter density around the Galactic center is 500 pc or less. 
This is a clear advantage of the ON/OFF method when 
compared to the background subtraction technique that 
is applied in JtJ which relies on the simultaneous obser¬ 
vation of the Galactic center region and a background 


region in the same finite H.E.S.S. field of view with ~ 2° 
radius. 

The application of standard quality criteria for H.E.S.S. 
data [Qj and the additional requirement for compati¬ 
ble instrumental and atmospheric conditions within an 
OFF1/ON/OFF2 observation result in a total of six 
OFF1/ON/OFF2 datasets. All datasets were taken 
within one week in 2010 with the H.E.S.S. I array of 
four identical IACTs. The total dead-tinre corrected ob¬ 
servation time for each of the three observed regions is 
3.05 h. The mean zenith angle of the array pointing for 
the datasets is 12 °. 

DATA ANALYSIS 

The image cleaning (see i) low and high pixel inten¬ 
sity thresholds for the data are chosen to be 7 pe (photo 
electrons) and 10 pe. Using the observed distribution 
of pixel intensities in cosmic ray events, it was checked 
that these image cleaning cut criteria eliminate effects 
due to differences in sky brightness between the observed 
regions. Standard Hillas criteria [§| for the selection of 7 - 
ray events are applied to the data. The thresholds used 
for image cleaning lead to an energy threshold of 290 
GeV. Only events with reconstructed directions within 
the central 2 ° angular distance around the pointing posi¬ 
tion of each observation are considered to account for 
the truncation of 7 -ray images near the edges of the 
H.E.S.S. field of view. The Galactic plane (| 6 | < 0.3°) 
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is excluded from the analysis to avoid the detection of 
7 -rays from astrophysical sources (e.g. the Galactic cen¬ 
ter source HESSJ1745-290, 18|) without relation to dark 
matter annihilation. The exclusion region is shifted by 
the respective pointing position offset in right ascension 
into the two background regions to equalize the accep¬ 
tance in the signal and background regions (see Fig. 0 . 
To rule out the detection of 7 -rays from astrophysical 
sources, the considered data with the chosen exclusion re¬ 
gions are analyzed with the ring background [l3 | method 
and a correlation radius of 0.1 deg prior to the ON/OFF 
analysis. The resulting skymaps of the three observed 
regions show no indication for a significant excess. It is 
concluded from the analysis with the ring background 
method that the chosen exclusion regions are sufficient 
to exclude astrophysical sources of gamma rays for the 
ON/OFF analysis. 

The mean exposure ratio, a = 0.5, for the ON/OFF data 
analysis is the ratio of the livetimes for the observation of 
the signal and background regions 0 - However, imbal¬ 
ances in the acceptance for background events between 
the signal region and the two background regions lead to 
a systematic error, cr Q , on the exposure ratio. A conser¬ 
vative estimate for the relative systematic error on the 
exposure ratio, <r a /a = 2%, is derived. This estimate 
results from a comparison of the number of events which 
pass 7 -ray event selection criteria in the two background 
regions. 


RESULTS 


A total of Aon = 24268 signal and Noff = 49028 
background events are measured that pass standard 
Hillas criteria [§j for the selection of 7 -ray events. The 
total 7 -ray signal s has a statistical significance of — 0.5er. 
The statistical significance is calculated with the log- 
likelihood ratio test statistic as described in [ 10 ] with 
the likelihood function (see also 0 ) 

L = P(Nqn, ab + s)P(Nqff, b)G(d , a, a a ) . ( 1 ) 


Here, P and G represent the Poisson and Gaussian distri¬ 
butions. The parameters b (mean number of background 
events) and a (exposure ratio with mean a) are treated 
as nuisance parameters. For comparison, the significance 
of the 7 -ray event excess as calculated with Eq. (17) 


11 ] without consideration of the systematic error on 


the exposure ratio is — 1.3<r. Since no significant 7 -ray 
signal is measured, an upper limit on the integrated 7 - 
ray signal for energies ranging from the instrumental en¬ 
ergy threshold to a maximum energy E is derived. For 
the calculation of the upper limit, the likelihood function 
that is given by eq. 0 is analyzed with the method de¬ 
scribed in 0.' The upper limit on the energy integrated 
signal translates (see e.g. i) into an upper limit on 
the velocity averaged dark matter self annihilation cross 



FIG. 4. Upper limits on the velocity averaged dark matter self 
annihilation cross section as a function of the dark matter par¬ 
ticle mass. The upper limits for the cored Einasto and NFW 
density profiles hold for a core radius of 500 pc and the anni¬ 
hilation of dark matter particles into light quarks (i)- The 
filled area around the upper limit curve for the cored Einasto 
dark matter profile shows the ±1 <t variations around the up¬ 
per limit that is expected for this dark matter density profile 
when no annihilation signal is detected. The derived upper 
limit is stronger than the expected upper limit due to the 
negative significance of the measured excess. For comparison, 
the velocity averaged annihilation cross section of a thermal 
relic dark matter particle is shown. Additionally shown are 
the upper limits that are derived in [j] for cusped Einasto 
and NFW profiles as well as the upper limit that is derived 
in Q for a cored dark matter density distribution around the 
Sculptor dwarf galaxy. 


section, (av)(M), for a dark matter particle with mass 
M = E. The variation of the instrumental response with 
the zenith and azimuth angles of the array pointing and 
within the field of view is accounted for in the analysis. 
The consideration of the 2% relative systematic error on 
the exposure ratio increases the upper limit on ( av ) by 
a factor of ~ 3. Upper limits on (av) are presented in 
Fig. [I] for Einasto and NFW dark matter density profiles 
with a 500 pc radius core of constant dark matter den¬ 
sity around the Galactic center. The parameters for the 
NFW and Einasto density profiles are taken from fl2| . 
The derived upper limits on (av) hold for the 7 -ray en¬ 
ergy spectrum that is expected from the self annihilation 
of dark matter particles into hght quarks (see 0 , the 
same spectrum is assumed in |7j). For an Einasto dark 
matter profile that is cored in the inner 500 pc around the 
Galactic center, values of (av) ~ 3 • 10 - 24 cm 3 /s or larger 
are excluded for dark matter particle masses in between 
~ 1 to ~ 4 TeV at 95% CL. The upper limits on (av) that 
are derived for an Einasto dark matter density distribu¬ 
tion with a core radius of 500 pc are the most constrain¬ 
ing exclusions that are derived for TeV mass dark matter 
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R (kpc) 

TEinasto 

^ON 

TEinasto 

'-'OFF 

tNFW 

J ON 

tNFW 

JOFF 

0 

2167 

268 

559 

78 

0.5 

1036 

268 

256 

78 

0.75 

636 

268 

165 

78 

1 

426 

255 

117 

75 

2 

138 

126 

46 

43 


TABLE I. Field of view averaged astrophysical factors for the 
signal (subscript ON) and for the livetime weighted average 
of the two background regions (subscript OFF). The values 
are in units of GeV 2 cm -6 kpc and are tabled for Einasto and 
NFW profiles as a function of the radius ( R ) of the central 
dark matter core. 

without the assumption of a centrally cusped dark matter 
density distribution in the search region. However, these 
limits are one order of magnitude less constraining than 
the current best limits for cusped dark matter density 
distributions (see Fig. 11 and two order of magnitudes 
weaker than the expectation for thermal relic dark mat¬ 
ter (see e.g. 0 )- 

For core radii different from 500 pc, the upper limit on the 
velocity averaged dark matter self annihilation cross sec¬ 
tion scales like (crv) r = (AJ 50 o P c/AJr) (av) 5 oo P c where 
A J denotes the difference between the held of view aver¬ 
aged astrophysical factors in the signal and background 
region and the subscript is equal to the core radius. The 
field of view averaged astrophysical factors in the signal 
and background region of the considered ON/OFF anal¬ 
ysis for different core radii are listed in table |T| The upper 
limits on ( av ) increase by a factor of 2 (5) if the radius 
of the central core of constant dark matter density is 750 
pc (1 kpc) when compared to a core radius of 500 pc. 


SUMMARY 

A search for a signal from annihilating dark matter 
around the Galactic center was performed. For this pur¬ 
pose, data that were acquired in dedicated ON/OFF ob¬ 
servations of the Galactic center region with H.E.S.S. 
were analyzed. No significant signal was found. The 
employed observation technique enabled the derivation 
of upper limits on (av) that are significantly more con¬ 
servative in respect to the distribution of dark matter in 
the Galactic center region than previous constraints. In 
particular, the constraints apply also under the assump¬ 
tion of a core of constant dark matter density around the 
Galactic center. If the dark matter density in the central 
500 pc around the Galactic center is constant and fol¬ 
lows outside of the core radius an Einasto profile, values 


of (av) that are larger than 3 • 10” 24 cm 3 /s were excluded 
for dark matter particle masses between ~ 1 and ~ 4 
TeV at 95% CL. This is currently the best constraint on 
(av) that has been derived without the assumption of a 
centrally cusped dark matter density distribution in the 
search region. 
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